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Abstract. In this article, we continue with our investigation of the Diophan- 
tine equation ~ = i + ■= and in particular its number of solutions R(n; a) for 
fixed a. We prove a couple of mean value theorems for the second moment 
(R(n;a)) 2 and from which we deduce log R(n; a) satisfies a certain Gaussian 
distribution with mean log 3 log log n and variance (log3) 2 log log n, which is 
an analog of the classical theorem of Erdos and Kac. And finally these results 
in all suggest that the behavior of R(n; a) resembles the divisor function d(n' 2 ) 
in various aspects. 



1. Introduction 

In the previous memoir of this series (see |HVj ) we studied the mean value 

S(N;a)= R ^ a ^ C 1 ) 

n<N 
(n,a) — 1 

of the number a) of positive integer solutions to the Diophantine equation 

all 

- = - + -■ (2) 

n x y 

Here we extend our investigation to the second moment and some consequences 
thereof. 

Theorem 1. For fixed positive integer a, we have, for every N G N with N > 2, 

2 

R(n; a) — 'TTT *(-») 

X mod a u\n 2 



E 

n<N 
(n,a) — l 



« a iVlog 2 TV, 



where <^ a indicates that the implicit constant depends at most on a, and where \ 
denotes the principal character modulo a. 

In the character sum here the term x = Xq contributes an amount d(n 2 ) where d 
is the divisor function and we can expect that this is the dominant contribution on 
average. Thus as a consequence of the Erdos-Kac theorem, just as for the divisor 
function d(n), one can anticipate that logi?(n;a) has a Gaussian distribution. As 
a first approximation we establish the normal order of log R(n; a). 
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Theorem 2. When a is fixed, the normal order of log R(n; a) as a function of n 
is (log 3) log log n. 

Let 

1 f z t 2 
$0) := —= I er^dt. 

V27T J-oo 

Then with a little more work we can establish the full distribution. 

Theorem 3. For fixed positive integer a, we have 

1 , f ^ T log R(n; a) — (log 3) log log n 1 , , . 
lim —card < n < N : & \ ' \ \ / — < z } = $ z . 
n^oo N \ ~ (log3)Vloglogn - J 

For completeness we also establish the mean square of R(n; a) for fixed a. Since 
R(n; a) resembles quite closely the divisor function d(n 2 ) in many aspects, we expect 
that their mean squares share the same order of magnitude. Thus the following 
theorem can be compared with the asymptotic formula 

d 2 {n 2 ) = N P & {\ogN) + OiN 1 - 5 ) 

n<N 

which holds for some 5 > and with Ps(-) a polynomial of degree 8. 
Theorem 4. Let a be a fixed positive integer and e > 0. Then 

R{n;a) 2 = NP 3 (logN;a) + O q (7V 35/54+£ ) 

n<JV 
(n,a) — 1 

where Pg(- ; a) is a degree 8 polynomial with coefficients depending on a, and its 
leading coefficient is 

6 1 \ / l x " 



s^nf'-iVni'-- ->"-- 



8!a 2 \ p J \ P P 2 J \ P 

p\a p\a 

The error term in the theorem above is closely related to the generalised divisor 
problem, and in particular depends on a mean value estimate for the ninth moment 
of Dirichlet L-functions L(s,x) modulo a inside the critical strip. As is easily 
verified, the error can be improved to O a {N 1 / 2+e ) under the assumption of the 
generalised Lindelof Hypothesis. 

2. Proof of Theorem [T] 

We rewrite equation ([2]) in the form 

(ax — n)(ay — n) = n 2 . 

After the change of variables u — ax — n and v — ay — n, it follows that R(n; a) 
is the number of ordered pairs of natural numbers u, v such that uv — n 2 and 
u = v = — n (mod o). 

Under the assumption that (n, a) = 1, R(n; a) can be reduced further to counting 
the number of divisors u of n 2 with u = —n (mod a). Now the residue class 
u = —n (mod a) is readily isolated via the orthogonality of the Dirichlet characters 
X modulo a. Thus we have 

R ( n '^ = 6(a) ^ X{~n)^2x(u), (3) 

^ ' X u \n 2 
mod a 
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where the condition (n, a) = 1 is taken care of by the character xi n )- 
Hence 

2 



E 

n<N 
(a,n) = l 



X mod a u\n 2 



X =X 







/AT 



n=l 
(a,n) — 1 



E ^nE^m 



X mod a 



xVx 



u\n 2 







= E E X 1 X 2 (-l)EEE^i(^2WXiX 2 (n)e-/ JV , 

X^mod a X2 mo< ^ a 1 u|n 2 f|n 2 

X-^ 7^Xq X2 "?~Xq 

where x denotes the principal character modulo a. In order to evaluate the sum 
over n, we analyze the Dirichlet series 

00 

/x r x 2 ( s ) : = E EE^w^^i^w^ 8 - 

71=1 u\n 2 v\n 2 

The condition u\n 2 can be written as uiu||n 2 with u\ squarefree, i.e. UiM 2 |n, and 
likewise for v\n 2 . Thus 



f ( s ) - V X^M ^Fjd) 



m — 1 



(4) 



where 



F ( d ) = E M 2 («i)/« 2 («i)Xi(«iW2)x 2 (wiv|)XiX 2 (rf)- 



Ul,U2|VliV2 

d=[uiu 2 ,fif2] 



The function F is multiplicative and so the inner sum above is 

n i+ef(?v s , 



k=i 



where 



F (p k )= E A* 2 («i)M 2 ("i)Xi(«i^)X2(«i"2)XiX 2 (P fc )- 

[lilt 

In particular we have 



(5) 
(6) 



Ul,U2,fl,f2 
[MlM2,fl'U2]=p fc 



^(p) = x (p) + E 

xeA-\{ Xl x 2 } 

where A" = {x x , X 2 > Xi X 2 > Xi X 2 , X x , X 2 , X x X 2 , X x X 2 } ( and the entries are considered 
to be formally distinct), and 

\F(p k )\<8k. 

Thus the Dirichlet series / converges absolutely for a > 1 and 

fx 1 ,x 2 (s) = G Xl , X2 (s)L(s, Xo ) I] L M, (7) 



4 



JING-JING HUANG AND ROBERT C. VAUGHAN 



where G x ^ (s) is a function which is analytic in the region -fts > 1/2 and satisfies 

G(s) < 1 (a>± + 6) 

for any fixed 6 > 0. As X : ,X2 are n °t characters of order 1 or 2, /x r x 2 ( s ) has a 
triple pole at s = 1 when x 1 = X2 or X1X2 = Xqi an d a simple pole otherwise. By 
Corollary 1.17 and Lemma 10.15 of |MVj . for fixed a, 

L(s ' x) -^iy <<2+|t| 

uniformly for a > ^ where -E(x) is 1 when x = x and otherwise. Hence by (5.25) 

of [MY] 

00 -. i-9+ioo 

EEExi(«)x 2 fe 2 We- n/N = ^- / Xi , X2 ( s )^r( s )d s , 



where 6 > 1. Since the gamma function decays exponentially fast on any vertical 
line we may move the vertical path to the f-line picking up the residue of the 
integrand at s = 1. The residue contributes an amount 

< iV(logA0 2 

and the new vertical path contributes 

This completes the proof of Theorem [T] 

3. Proof of Theorem [2] 
By Theorem [U we expect that for almost all n with (a, n) = 1, R(n; a) is close to 

E X(-n) £*(«)• 

X niod a it | n 2 

X =Xq 

Thus we need to examine the contribution from the characters modulo a of order 
1 and 2. For general a, there may be many quadratic characters modulo a. Nev- 
ertheless we believe that the major contribution to the sum above comes from the 
principle character, and this is of size 

d{n 2 ) 
0(o) ' 

Thus, for fixed a, log R(n; a) should have the normal order of logd(n 2 ), namely 
(log 3) log log n. When (n, a) > 1 we have 

R(n; a) — R[n/ (n, a); a/ (n, a)) (8) 

and so we can expect that the general case follows from the special case (n, a) = 1. 

Before embarking on the proof of Theorem [21 we state a lemma. We define, for 
any quadratic character x, 

fl x (n) = card {p, k : k > l,/|n,x(/) = l} . 
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Lemma 5. Suppose that \ is a quadratic character to a fixed modulus a and that 
N > 3. Then 

J2 (r> x (n)-iloglogAM <jVloglogiV 

n<N \ ~ J 

and 

(fix(n)-^loglogn) <iVloglog7V. 

Kn<N \ / 

Proof. The proof follows in the same way as Turan's theorem (see Theorem 2.12 of 
|MVj ) on observing that 

]T - = l\og\o g N + 0(l) 
* — ' V 2 

x(p)=i 

and this is readily deduced from Corollary 11.18 of |MVj . □ 

It is an immediate consequence of the above lemma that VL X (n) has normal order 
^ log log n. In particular, for any fixed e > 0, for almost all n, 

gf2 x (n) ^ 3(5+^) log log n 

Now, for any quadratic character \ modulo a, let 

9x( n ) = Y x( u )- 

u\n 2 

This is 

n {i+x(p)+x 2 ( P )+---+x 2k (p))- 

p k \\n 

When x(p) = — 1 the general factor is 1, and when x{p) = 1 it is 2k + 1. Hence 

< 9 x (n) < 3 n *H 

Thus for any fixed e > 0, for every quadratic character modulo a, for almost all n, 

g x {n) < (log»)(4 ]o « 3+e ). (9) 

Let 

r ( n;a ) = 7T77 x(-n/(n,a))g x (n/(n,a)). 

4>(a/(n,a)) 

^ W V ' X mod o/(n,o) 

x 2 =x 

Since R(n;a) — R(n/(n,a);a/(n,a)), it follows by Theorem Q] that 

^(i?(n;a)-r(n;a)) 2 ^ (i?(m; a/rf) - r(m; a/d)f < iV(log iV) 2 . 

n<JV <2|a m<N/d 

(m,a/d) — l 

Hence, for any fixed e > 0, for almost all n we have 

\R(n;a)-r(n;a)\ < (logn) 1+e . 
Therefore, by ([5]), for almost all n, 

d({n/{a,n)Y\ 



R(n; a) 



(f>(a/(a,n)) 



<(\ogn) L+2e . (10) 
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Now 3 < d(p 2k ) = l + 2k<3 k . Hence 

3 W (n)-*(a) < d ((n/( a , n)) 2 ) < 3° (n) (11) 

and it follows that 

, „ d ({nl{a, n)) 2 ) , Q 

(\ogn) l °^ < V ;, >> > < (logn) 1 ^^ 

for almost all n. Theorem [5] now follows. 

4. Proof of Theorem [3] 
By (ITO and (fTTj) . for every fixed e > 0, for almost all n, 

5u(ti) 

— - -, ^ - (logn) 1+£ < R{n;a) < + (log,i) 1+£ . 

(p(a/{a,n)) 

Moreover, for almost all n we have f2(n) > u>(n) > (1 — e) log log n. Hence for any 
8 with < 8 < log 3 — 1 we have, for almost all n 

3 u(n)-^(a)-lo g <p(a/(a,n)) ex ^_^ ogn yS^ < fl ) < 3 0(n) eX p((log n)- S ) 

and so 

3 u < n > exp(-eVloglog^) < a) < 3 n(n) exp(eVloglogn). 



Let 



and 



S(N; z) — card in < N : l0g *<" ; a) ~ (1 ° g 3) l0g l0g " < z 

log 3 V log log n 

.s. ! V::) = n ! r,l><iV: t]( " HOgl0g "<z 

V log log n 

V ,.V: :» = ,„, .Nn^iV:^- l0gl ° gn <4. 

Vl°g log n 



Then for a non- negative monotonic function tending to sufficiently slowly as 
N — > oo we have 

-T}(N)N + S-(N; z-e) < S(N; z) < n(N)N + S+(N; z + e). 

Hence, by the Erdos-Kac theorem (see, for example Theorem 7.21 and Exercise 
7.4.4 of [MY]), 

$(z-e) < liminf N- 1 S{N;z) < lim sup N~ 1 S(N; z) < ®(z + e). 
The theorem now follows from the continuity of $. 

5. Proof of Theorem [4] 

By a similar discussion to that in $2j we can show that the generating Dirichlet 
series for R(n; a) 2 is 

R{n;a) 2 1 ^ _ 
E —^ = J(a^ E X 1 X 2 (-l)f XvX2 (s), 

n=l rv ' Xi'*2 

(n,a) = l mod a 

where fx 1 ,x 2 i s ) ^ s analytic in the region 3?s > 1/2 and is given by (0. Here 
fx 1 ,x 2 ( s ) nas a P°l e a * 1 °f or d er at least 1, and as high as 9 exactly when \i and 
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X 2 are equal to the principle character x ■ Now on applying Perron's formula, we 
have for = 1 = 1 + 1/ log(2 JV) ; 

n<7V Xi 
(n,a) — l mod a 

(12) 

Since we are shooting for the asymptotics for the the mean square, smoothing 
factors of the kind used in section [5] are best avoided. Since the integrand includes 
a product of nine L-functions, we cannot expect to to be able to move the vertical 
integral path too close to the 1/2-line in the current state of knowledge. Neverthe- 
less, the following result of Meurman |Melj which extends Heath-Brown's theorem 
[H-B] on the twelfth power moment of the Riemann zeta function to Dirichlet L- 
functions, provides a starting point for the analysis. 

Lemma 6. 

rT 



£ / |L(± + zi,x)r*«a 3 T 2 ^, 



X mod a 

where e > 0, a > 1 and T > 2. 

Then, adapting the argument of Chapter 8 of Ivic [Iv] for the Riemann zeta 
function to Dirichlet L-functions establishes the following. 



Lemma 7. 



f \L(§ + it, x )\ 9 dt<: a T 1+s , 

J-T 



-T 

where e > 0, x is a fixed Dirichlet character modulo a > 1 and T > 2. 

If one utilizes the sharpest estimates for the underlying exponential sums, Lemma 
[7] is susceptible to slight improvements. 

Now, we move the vertical integral path in (TT2"j) to the 35/54-line, picking up the 
residue of the integrand at 1 . Thus 

e+iT N s /.35/54-iT N s r-e+iT N s 

fx-,,x 2 { s )— ds = fx r x 2 ( s )— ds + fx r x 2 ( s )— ds 

-iT 1 i s Je-iT 1 A s J35/54+iT 1 1 3 

1-35/54+iT N s / N s\ 

+ / fx, ,x 2 (s)—ds + Res s=1 f x , x (s) — 

Here, in order to deal with the contribution from the horizontal integrals, we 
cannot afford to use the crude convexity bounds on Dirichlet L-functions, due to 
the large number of L-functions in the integrand. Fortunately, a sharper bound 
has been established by Pan & Pan in Theorem 24.2.1 of |PP) . 



Lemma 8. Let I > 3, L = 2 l 1 and o\ = 1 — l(2L — 2) 1 . Then when a > o\ 

L(a + i t, X )«a \t\ x l^ L -^ log |t| 
holds uniformly for \t\ > 2. 
When / = 3 we obtain 

L(a + it, X )<Za \t\ 1/e log\t\ 
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uniformly for |i| > 2 and a > \, and when I = 4, 

L(a + it,x)<^a |t| 1/14 log |t| 
uniformly for a > 5/7. Thus, by the convexity principle for Dirichlet series, 

L(a + it :X )^a |<r (CT)+e 
uniformly for \t\ > 2 and a > ^ where 



I - !("-!) wheni<a<|, 
M<» = <( ^ when f < cr < 1, 

when 1 < a. 

We note that = J| < ± and M (f) = ^ 

Now the horizontal paths contribute 

« / ^ N°\f x x {a + iT)\T- l da. 

J35/54 1 1 

and this is 

< max N v T 9^)-i+e^ 

35/54<<T<l+e 

and by the piecewise linearity of a and this is 

<C N 1+e T~ 1 + 7V 5 / 7 j l9 ^( 5 / 7 )^ 1 + £ + ^35/54^9^(35/54) 

When T = N this is 

< jy35/54+<^ 

On the other hand, by Lemma [7] the vertical path also contributes 

< Ar 35/54+e_ 

The main term comes from the residual contributions, which, in the case that 
X 1 = X 2 = X > i s ATPg (log N; a) where P 8 (- ; a) is a polynomial of degree 8 whose 
coefficients depend on a. Notice that for other choices of \i an d X 2 ' ^ ne res idual 
contribution gives a polynomial of log N of lower degree than above. 

For the leading coefficient, we need more precise information about fx ,x - By 
Q, © and © we have 

/x ,x = Ha, x ) II f 1 + E = M«. X ) 9 + 6 ^ S + ^ 2S )(! - ^ S ) 6 > 

p\a \ k=l / pfa 

from which the leading coefficient is readily deduced. This completes the proof of 
Theorem |U 

In conclusion we remark that a concomitant argument will give 
J2 R(n; a) k = NP&^Qog N; a) + O a (N a ^) 

n<N 
(n,a) — 1 

for any e > 0, where P 3 k_ 1 (- ; a) is a polynomial of degree 3 k — 1 whose coefficients 
depend on a and a*; is a constant that depends on the best 3 fc -th power moment 
estimates for L(s, x) in the critical strip and the quantity n(cr) defined above. This 
question is closely related to the generalised divisor problem, and one is referred to 
Chapter 13 in Ivic [Iv] for more details. 
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